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Abstract 

We discuss recent developments in lattice QCD for the bulk thermodynamics of the transition. 
We review the current status of the equation of state, the possible implications of a magnetic field 
and the fluctuations of conserved charges, like the net electric charge or baryon number. We also 
show predictions for higher cumulants, which will be experimentally available in near future. 



1. Introduction 

Lattice simulations aiming to describe the features of quark gluon plasma and its transition to 
confined matter have become a field of high precision in the past years. An important goal of the 
efforts on the lattice is to provide theoretical background for collision experiments. The relevance 
of lattice results are most prominent for the heavy ion physics pursued at RHIC and LHC, where 
matter with high energy density, but moderate or small chemical potentials are produced. 

Lattice simulations solve Quantum Chromodynamics (QCD) in equilibrium. The solution 
emerges from a continuum limit, which is an extrapolation from a set of simulations on lattices 
with smaller and smaller lattice spacings (a). Finite temperature (T) is naturally implemented 
by a torus in Euclidean time. In most thermodynamics studies the lattice spacing is expressed 
through the number of lattice points in the Euclidean time direction: N,, which translates to 
a = l/(N,T). In most early papers N, was set to 4 or 6. By today N, = 16 data sets are not 
uncommon, and are believed to be necessary for a controlled continuum extrapolation even with 
improved actions. 

Longer correlations (inverse pion mass) require larger computer resources to simulate the 
theory. In the past decade lattice field theory has seen a great development in gradually reducing 
the quark masses down to their physical values [ 1 1. Simultaneously to this trend our increasing 
understanding of various discretization effects has lead to more efficient discretization schemes 
(actions) that allow a continuum extrapolation at an affordable cost 13 |3). In view of these 
developments here we only discuss results with physical (or nearly physical) quark masses close 
to or in the continuum limit. 

2. Transition temperature 

Evidence form lattice QCD shows that at zero chemical potential the transition from the 
chirally broken, confining phase of hadrons to the quark gluon plasma phase is a crossover (4). 
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The transition temperature of a crossover transition depends on the observable one chooses. 
At vanishing light quark masses the chiral condensate is a genuine order parameter, which can be 
used to extract a chiral transition temperature. At sufficiently high quark masses the deconfine- 
ment transition is identified by a discontinuity in the Polyakov loop as a function of temperature. 
Varying the quark mass in the space of theories the relevance of chiral and deconfinement aspects 
can be tuned. With sufficiently light u and d quarks an 0(4) scaling of the chiral observables is 
conjectured |5 , 6 1. This universal mass dependence has been used to relate the transition temper- 
ature (T c ) in QCD to T c for the actual simulation parameters Q. 

The Wuppertal-Budapest collaboration has simulated the QCD transition with physical quark 
masses in the stout staggered formulation [8 1. The lattice resolutions range from N, = 6 to N, = 
16. From the inflection point of the chiral condensate we find T c = 155(3)(3) MeV, where the first 
of the errors is statistical, the second is systematic, including the uncertainties of the scale setting 
OQTO. Other T c definitions from chiral observables span a range of 147 - 157 MeV. Recently 
the HotQCD collaboration has published a transition temperature, 154(8)(1) MeV, based on the 
0(4) scaling [7|, which is very well compatible with the range given by the Wuppertal-Budapest 
collaboration in 2006. This is in contrast to the earlier continuum result (192(4)(7) MeV) by the 
RBC-Bielefeld group. A key improvement since then was the adoption of the HISQ action Q. 

The quoted values of T c are all based on lattice observables that are remnants of the chiral 
order parameter. One can also study other observables, such as the Polyakov loop, which is 
the remnant order parameter of the deconfinement transition. The Polyakov loop as well as the 
inflection point is renormalization scheme dependent. E.g. the scheme used in Ref. was 
different than in our latest work ifTOl . We made this change in our convention to be compatible 
with Ref. IfTTI in every detail. While the comparison of our result to the asqtad data has showed 
a disagreement, the new HISQ data of the HotQCD collaboration was in agreement with ours |7|. 
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Figure 1: Left: The renormalized chiral condensate. We used the inflection point to define a transition temperature: 
155(3)(3) MeV. Right: the renormalized Polyakov loop. Although one can precisely measure its temperature dependence, 
one cannot easily pinpoint an inflection point. The shape of the Polyakov loop curve, as well as any derived transition 
temperature depends on the details of the renormalization scheme I 9] |10l . 

There are other observables that indicate deconfinement. The normalized energy density 
e/T 4 has an inflection point at 157(4)(3) MeV, for the trace anomaly we have 154(4)(3) MeV. 
The speed of sound has a minimum at around ~ 145(5) MeV. These temperatures cover roughly 
the same range as we find from chiral observables. From the strange susceptibility, which charac- 
terizes the deconfinement of the strange degree of freedom we have a somewhat higher inflection 
point. Of course, in a broad crossover the definition of the most singular point is also ambiguous. 

The concept of an identical chiral and deconfinement transition temperature was promoted by 
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Figure 2: Left: The trace anomaly at several lattice spacings, either with 2stout action of the Wuppertal-Budapest 
collaboration 1 1 5 16 1, or with the HISQ action of the HotQCD collaboration 1 17 1. For comparison we show the Hadron 
Gas Resonance model's prediction as well as the historic asqtad data of Ref. I'll]. The use of the HISQ action has 
significantly reduced the discrepancy between the two collaborations. Right: The speed of sound by the Wuppertal- 
Budapest collaboration at three lattice spacings 1 15] (these data are not yet available from HISQ calculations). 



early lattice data [ 12 1, that positioned the peak of the chiral and the Polyakov loop susceptibilities 
at the same temperature: T c , Lacking resources, only theories with heavier quarks were simulated 
at the time. Results with somewhat smaller quark masses on coarse lattices confirmed this picture 
iTPJl . Solutions to two- flavor QCD based on the functional renormalization group have shown 
the same conclusion Ifl4ll . In full (2+1 flavor) QCD, however, we see a very broad transition. 
Thus the transition temperatures based on various observables generally do not coincide, they 
have their values in this broad transition region. 



3. Equation of state 

The pressure and energy density of the quark gluon plasma is one of the key targets of lattice 
QCD simulations. These form the equation of state, which is required to close the equations of 
relativistic hydrodynamics IfTSll . 

Early works on the QCD equation of state were dominated by quenched results Ifl9l l20l . 
These were recently extended in range and precision, so that a connection to perturbation theory 
could be made [21 j. The calculations were also extended towards larger gauge groups [22 J. For 
several years we have seen various exploratory studies with non-physical or not renormalized 
quark masses. Since Ref. |8| physical simulations have been feasible. In staggered simulations 
discretization effects raise the mass of the hadrons, and through this T c appears to be higher, and 
the energy and entropy density function is found steeper than physical, and the trace anomaly 
shows a higher peak. This is an artefact ("taste-breaking") which is absent after continuum 
extrapolation, its impact on coarse lattices and older actions (like asqtad) was emphasized in 
Refs. [9 10 23). The Wuppertal-Budapest collaboration exploited the favourable taste-breaking 
features of the stout action and calculated a continuum estimate of the equation of state lfl5l 
based on N, - 6,% and 10 runs. HotQCD used the asqtad and p4 actions at N t = 6 and 8 
ifTTl [24), though these were later found to be incompatible with results with the better HISQ 
action [ 17 1. For a comprehensive review see Ref. J25). 

At this conference, we heard two contributions discussing the equation of state. The Wuppertal- 
Budapest collaboration has updated its continuum estimate of 2010 [ 15 1 by adding a finer dataset 
(N, = 12) and an even finer check point (N t = 16) and using these in a continuum extrapolation 
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Figure 3: The entropy density as a function of temperature. Left: We show the lattice data by the Wuppertal-Budapest 
collaboration (N, = 6, 8 and 10) together with the parametrization (solid red line) of these data sets published in Ref. 1 15]. 
The black dots show the result of a recent continuum extrapolation based on increased statistics and finer lattices (N, = 
12). We see almost no deviation from the earlier parametrization. Since there is no HISQ data for the entropy density yet, 
we plotted HotQCD's asqtad data set at N, = 8 for comparison (magenta triangles) II II . As discussed already (Fig. [5), 
this set features a higher-than-physical transition temperature and thus it shows steeper equation of state and higher peak 
in the trace anomaly. As remedy, the s95p parametrization was introduced 1 23 1 by replacing lattice data in and below the 
transition range by the Hadron Resonance Gas prediction. By this, the slope of the entropy curve was left unchanged, but 
the whole curve was shifted up towards the Stefan-Boltzmann limit. Right: The contribution of the charm sea quark to 
the entropy density. The Wuppertal-Budapest collaboration's data with N t = 6, 8 and 10 is shown for dynamical 2+1 + 1 
flavor simulations together with its charmless continuum data and the respective perturbative estimates 1161 . 



Ifl6l . A. Bazavov (HotQCD collaboration) has presented data at N, = 4,6,8, 10 and 12 IfTTl . 
First, the trace anomaly was obtained form the lattice simulations. Then, the other quantities, 
like pressure, entropy density or speed of sound can be derived using thermodynamic relations. 
The results are summarized in Fig. [2] and [3] The missing contribution of the charm quark has 
been estimated in C. Ratti's talk, by actually performing new 2+1 + 1 flavor simulations (see 
Fig. fright). The charm's contribution to the pressure or energy density was found to be relevant 
from about 250-300 MeV temperature ifTBI . 

4. Phase diagram 

The transition line in the T - /j.b phase diagram at small chemical potentials is necessarily a 
line of crossover, with different width and curvature for various observables. Direct Monte-Carlo 
simulations of the four-dimensional quantum field theory at finite baryo-chemical potential are 
at present not possible. Nevertheless, starting with Ref. [29 1 a renewed interest has been seen for 
fiB > questions in lattice QCD, and several indirect strategies exist BUI [3D . 

The curvature of the transition line in full QCD has been found very small by recent lattice 
simulations. In Ref. [6| simulations on N, = 8 lattices were used to map the chiral condensate 
and its fx 2 B dependence onto a universal 0(4) scaling behaviour as a function of the quark mass. 
Evaluating the fits in the physical point gave k = 0.0066(2)(4), with k = -T c dT c (fi 2 B )/diJ 2 B \ =Q . 
The continuum result has been calculated by the Wuppertal-Budapest collaboration |26|, with 
simulations explicitly in the physical point. Their result was k = 0.0066(20). In this work, a 
similar analysis have been made for the strange susceptibility, there the curvature was found to 
be 0.0089(14). In addition, the width of the transition has also been extrapolated to hb > 0. This 
showed that the transition is not getting stronger with increasing chemical potential, but rather it 
stays a crossover (see Fig.^eft). 
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Figure 4: Left: The QCD phase diagram at small chemical potentials [26]. The transition temperature is plotted to 
leading order in Ty was extracted from the chiral condensate, and T*' is derived from the strange susceptibility. 
The solid bands represent the statistical uncertainty, while the color bands stand for the measured width of the transition. 
At this level of the expansion, there is no indication for the transition getting stronger with increasing fig. The dashed 
line show the chemical freeze-out curve of Ref. |27|. Right: Like T c , the equation of state can also be extrapolated to 
finite /ifi. Here /i^ refers to the baryo-chemical potential while the mean strageness was kept zero to linear order. For a 
detailed discussion of these continuum results and a simple parametrization, see Ref. 1 28 ]. 



This latter conclusion is in line with the work by Philipsen and de Forcrand claiming the 
absence of a critical endpoint, at least in the range where analytical continuation or Taylor ex- 
pansion can work l32l . This result was often confronted with the critical point found by Fodor 
and Katz 11331 . Note, however, that at small chemical potentials even Ref. [ 33 1 finds a slightly 
weakening transition, this trend then reverses at higher where the available Taylor expansion 
from imaginary /is no longer works. As the authors also emphasize, these results 021 33) came 
from coarse lattices. A continuum limit has not yet been feasible (because of the high computa- 
tional costs), which leaves us without any conclusive lattice evidence for a critical endpoint or a 
first order line at high yt/g. 

Having no singularity in a large range for chemical potentials, the equation of state, too, can 
be extrapolated to fi B > 0. By integrating the trace anomaly we have logZ at the \i B = 0. The 
baryon number susceptibility ~ <9 2 logZ/(5/VB) 2 | /iB= o is also accessible from lattice, thus, a linear 
expansion is possible. The continuum result was given in Ref. Il28l (see Fig.|4]i. 

We finally discuss the effect of a magnetic field. The strong magnetic fields in non-central 
heavy ion collisions [ 36 1 motivated the study of the QCD transition in the presence of an elec- 
tromagnetic background. A B field is easily implemented on the lattice 11371 [38l . Based on 
models (e.g. PNJL (39)) and exploratory lattice simulations [37] the expectation was to see an 
increase in T c with growing magnetic field. As a joint effort of the Budapest, Regensburg and 
Wuppertal lattice groups T c was calculated with physical quark masses in the continuum limit. 
T c was shown to decrease with B (see Fig.^left), while the transition was still a crossover 11341 . 
To show the role played by the quark mass, an auxiliary data set was taken with higher hadron 
masses. Then the trend was reversed, and the earlier (apparently contradicting) lattice result 071 
was reproduced. To explain the behaviour seen in the PNJL model the same group calculated the 
chiral condensate at B > and T > 0, and highlighted the differences between lattice data and 
the model's assumptions 0311 (see Fig. fright). 
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Figure 5: Left: The inflection point of the chiral condensate and of the strange susceptibility. The cross-over temperature 
decreases with a magnetic field, independently of its definition 1 34]. Right: The behaviour of the chiral condensate as a 
function of the magnetic field at various temperatures |35|. The magnetic catalysis at low T turns into a non-monotonic 
behaviour in transition region. 
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5. Fluctuations of conserved charges 

Correlations and fluctuations of conserved charges have been proposed long ago to signal 
the transition 14011411 . At high temperatures fluctuations are expected to approach the ideal gas 
limit. On the other hand, in the low-temperature phase they are expected to be small since quarks 
are confined and the only states with non-zero quark number have large masses. Agreement with 
the Hadron Resonance Gas (HRG) model predictions is expected in this phase. 

In the past year both the Wuppertal-Budapest [42] and HotQCD [43] collaborations published 
results for the second order fluctuations, i.e. the second derivatives of the free energy with 
respect to chemical potentials of various conserved charges, like baryon number, electric charge 
or strangeness. 

The lattice results are always normalized to volume and the respective power of temperature. 
E.g. the diagonal charge and the off-diagonal charge-baryon correlators are defined by 

Q= 1 d 2 logZ BQ= 1 d 2 logZ 

Xl vr 3 (dn Q m 2 ' Xu vr 3 (dfx a /T)(dfx B /T) ' 

The experimentally more interesting observables are ratios of these derivatives, where the 
unknown volume factor cancels. To describe the freeze-out parameters (T and /ub) specific ratios 
have been introduced and calculated (e.g. skewness / mean of the net charge yield) using the 
coarser data sets of the HISQ action (N, = 6, 8) in Ref. l43l . 

The charge kurtosis x variance could also be used as a thermometer, and was in the focus of 
C. Schmidt's talk |45l . In Fig.|7]we plot the temperature dependence ofx®IX® at zero chemical 
potential from both the BNL-Bielefeld group and the Wuppertal-Budapest collaboration. While 
we witness steady progress in the analysis of experimental data towards the precise measurement 
of second and higher order fluctuations [46, 47 48 1, lattice methods are increasingly successful 
in higher derivatives, which will finally enable the extrapolation of the kurtosis data towards 
finite chemical potential. 

6. Summary 

We discussed a selection of the recent developments in bulk thermodynamics, focusing on 
those results with physical quark masses close to or in the continuum limit. For technical reasons, 
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Figure 6: Diagonal and off-diagonal correlators of conserved charges as a function of temperature. Only continuum 
extrapolated lattice data are shown. There is a good agreement between HotQCD 1431 and Wuppertal-Budapest 1421 
data, obtained using the HISQ and 2stout actions, respectively. The broader temperature region by the Wuppertal- 
Budapest collaboration allows a comparison with the HRG model's prediction, and was used to benchmark resummation 
schemes at high temperature 1441 . 



all these results were obtained using the staggered formalism. The independence of the results 
with respect to the applied discretization scheme is an important consistency requirement of 
lattice QCD, and to this end there have been recent checks with overlap l49l . domain-wall ll50l 
and Wilson [51 1 fermions, with the compromise of using heavy pions in these comparisons. 

By today we know several aspects of the QCD transition through Euclidean-time observables 
at zero and small chemical potentials. Now the challenge for lattice QCD is to develop new tech- 
niques for extracting real-time physics (e.g. transport coefficients) and to increase the accessible 
range in fig, and eventually, to map the phase diagram. 



References 

[1] C. Hoelbling, PoS LATTICE 2010 (2010) 01 1 

[2] C. Morningstar and M. J. Peardon, Phys. Rev. D 69 (2004) 054501 

[3] E. Follana et al. [HPQCD and UKQCD Collaborations], Phys. Rev. D 75 (2007) 054502 

[4] Y. Aoki, G. Endrodi, Z. Fodor, S. D. Katz and K. K. Szabo, Nature 443 (2006) 675 

[5] S. Ejiri, et a/.Phys. Rev. D 80 (2009) 094505 

[6] O. Kaczmarek, et a/,Phys. Rev. D 83 (201 1) 014504 

[7] A. Bazavov, et a/.Phys. Rev. D 85 (2012) 054503 

[8] Y. Aoki, Z. Fodor, S. D. Katz and K. K. Szabo, JHEP 0601 (2006) 089 

[9] Y. Aoki, Z. Fodor, S. D. Katz and K. K. Szabo, Phys. Lett. B 643 (2006) 46 Y. Aoki, S. Borsanyi, S. Durr, Z. Fodor, 
S. D. Katz, S. Krieg and K. K. Szabo, JHEP 0906 (2009) 088 

[10] S. Borsanyi et al. [Wuppertal-Budapest Collaboration], JHEP 1009 (2010) 073 

[1 1] A. Bazavov et al, Phys. Rev. D 80 (2009) 014504 

[12] For an early review of lattice thermodynamics see: F. Karsch, Lect. Notes Phys. 583 (2002) 209 |hep-lat/01 06019 1. 

[13] M. Cheng, era/.,Phys. Rev. D 77 (2008) 014511 

[14] J. Braun, L. M. Haas, F. Marhauser and J. M. Pawlowski, Phys. Rev. Lett. 106 (201 1) 022002 

[15] S. Borsanyi et al. JHEP 1011 (2010) 077 

[16] C. Ratti (Wuppertal-Budapest collaboration). Proceedings of Quark Matter 2012, Washington D.C. 

[17] A. Bazavov (HotQCD collaboration), Proceedings of Quark Matter 2012, Washington D.C. |arXiv:1210.6312| 

[18] P. F. Kolb and U. W. Heinz, In *Hwa, R.C. (ed.) et al.: Quark gluon plasma* 634-714 

[19] G. Boyd et al. Nucl. Phys. B 469 (1996) 419 

[20] M. Okamoto et al. [CP-PACS Collaboration], Phys. Rev. D 60 (1999) 094510 

[21] S. .Borsanyi, G. Endrodi, Z. Fodo r, S. D. Katz and K . K. Szabo, JHEP 1207 (2012) 056 

[22] See e.g. B. Lucini and M. Panero, arXiv:1210.4997 [hep-th], and references therein 

7 



HRG 

BNL-Bielefeld preliminary N.=6 ^ • -h 
N t =8 „ j- . 

Wuppertal-Budapest preliminary N.=6 i — e — . 

N.=8 



N,=10 
N,=12 
N. 16 



140 160 180 
T [MeV] 



T d 4 log Z 

0.8 - 



0.4 



charge 1 — e — i 
baryon > — a — i 
light quark > — v — ■ 
strange quark 1 — e — 1 

Wuppertal-Budapest 
preliminary continuum 



0.2 

q 1 ^AAAAAAAAAAAAAAA A A A A A A A A A A AAAA AAAAAAAAAA ^ 

200 220 240 260 280 300 320 340 360 380 400 
T [MeV] 



Figure 7: Left: The normalized electric charge kurtosis from lattice QCD: the BNL-Bielefeld collaboration has used 
the HISQ action on N, = 6 and 8 lattices 1451 . the Wuppertal-Budapest collaboration shows results at N, = 6, 8, 10, 12 
and 16. From the plot we see that the use of fine lattice is essential for a controlled continuum limit. At present data 
indicates deviations from the Hadron Resonance Gas model's prediction even below T c . Right: At high temperatures 
the Wuppertal-Budapest collaboration has attempted a preliminary continuum extrapolation for xa for electric charge, 
strangeness and baryon number. The difference we see between light and strange quark's behaviour indicates that the 
mass of the strange quark is still noticeable up to about 260 MeV temperature. 
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